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ABSTRACT 



Context. In recent years, several potentially habitable, probably terrestrial exoplanets and exoplanet candidates have been discovered. 
The amount of CO2 in their atmosphere is of great importance for surface conditions and habitability. In the absence of detailed 
information on the geochemistry of the planet, this amount could be considered as a free parameter. 

Aims. Up to now, CO2 partial pressures for terrestrial planets have been obtained assuming an available volatile reservoir and out- 
gassing scenarios. This study aims at calculating the allowed maximum CO2 pressure at the surface of terrestrial exoplanets orbiting 
near the outer boundary of the habitable zone by coupling the radiative effects of the CO2 and its condensation at the surface. These 
constraints might limit the permitted amount of atmospheric C0 2 , independent of the planetary reservoir. 

Methods. A ID radiative-convective cloud-free atmospheric model was used to calculate surface conditions for hypothetical terrestrial 
exoplanets. C0 2 partial pressures are fixed according to surface temperature and vapor pressure curve. Considered scenarios cover a 
wide range of parameters, such as gravity, central star type and orbital distance, atmospheric N 2 content and surface albedo. 
Results. Results show that for planets in the habitable zone around K-, G-, and F-type stars the allowed CO2 pressure is limited by 
the vapor pressure curve and not by the planetary reservoir. The maximum C0 2 pressure lies below the CO2 vapor pressure at the 
critical point of p alt =73.8 bar. For M-type stars, due to the stellar spectrum being shifted to the near-IR, CO2 pressures above p cl j t 
are possible for almost all scenarios considered across the habitable zone. This implies that determining CO2 partial pressures for 
terrestrial planets by using only geological models is probably too simplified and might over-estimate atmospheric CO2 towards the 
outer edge of the habitable zone. 

Key words. Planets and satellites: atmospheres, Planets and satellites: composition 



1. Introduction 

Given the difficulties and challenges of detecting sub-surface life 
on Earth, any life to be first discovered beyond our own solar 
system will most likely be restricted to the planetary surface and 
atmosphe re. This is the basis of the concept of the hab itable zone 
(HZ, e.g.. lDoldl964lHar3ll978irKasting et alJl"99l . The HZ is 
defined as the region around a star where a rocky planet with a 
suitable atmosphere can host liquid water on its surface, a con- 
dition motivated by the fact that all life as we know it requires 
liquid water. 

Several studies have implied that small, potentially rocky 
planets are common (e.g.jHoward et al.ll2010j 1 Wittenmver et all 
201lllBorucki et alj201 WMavor et alj201 lllCassan et al.ll2012l 



lHaghighipour et all 120 1 PL [Tinnev et~a"fll201 ll) which raises the 
possibility of habitable satellites around these planets. 

A simple criterion for the potential habitability of a planet, 
which is immediately accessible from the discovery data, is its 



equilibrium temperature, 
calculated by 



T = 

1 eq — 



(l-A)F 

4cr 



0.25 



T eq . The equilibrium temperature is 



(1) 



Gaidos et alJ 120121) . Hence, it is not unreasonable to assume 
that planets in the HZ of their central stars may also be rela- 
tively common. Indeed, some potentially habitable (candidate) 



have already been discovered (lUdry et al 


II2007L MavoretalJ 


20091 Borucki et alJ 201 ll Pepe et all 201 1 


iBonfils etalJl201lL 


Anglada-Escude et al. 12012. Delfosse et al. 


2012. Borucki etal J 


2012|). Also. Neptune- or Jupiter- like planets have been dis- 
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where A is the planetary albedo, F the stellar flux at the 
orbital distance of the planet and c r the Stefan - Boltz mann 
constant. As was discussed b y, e.g., ISelsis et all d2007l) and 
iKaltenegger & Sasselovl d201 ll) . a habitable planet should have 
Teq < 270 K to avoid a runaway heating of the surface and cor- 
responding loss of the complete surface water reservoir. For low 
values of r eq near the outer edge of the HZ (e.g., model calcu la- 
tions forGL 581 d suggest T eq ~ 190K, lvon Paris et al.ll2010l) , a 
massive greenhouse effect must be provided by the atmosphere 
to obtain habitable surface conditions. 

H2O is the most obvious candidate of radiatively active gases 
which could provide the necessary greenhouse warming. It pro- 
vides the bulk of the greenhouse effect on Earth. Furthermore, 
H2O is by definition present on the surface of a habitable planet. 
The H2O partial pressure in an atmosphere of a potentially hab- 
itable planet is controlled by evaporation (or sublimation) from 
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the surface reservoir, taking into account the water vapor pres- 
sure curve. Besides H2O, CO2 is usually considered the most im- 
portant greenhouse gas for the determin ation of the outer bound- 
ary of the HZ (e.g.. feasting et alJll993h . 

On Earth, CO2 is controlled by processes such as volcanic 
outgassing or rock weathering. To estimate CO2 partial pres- 
sures for terrestrial e xoplanets, up to now only geological mod- 
els were used (e.g.. lEMns-Tanton & Seagerl 120081 [Kite et al] 
120091 iKiteetai] 1201 lL lEdson et al.1 120 12L I Abbot et all l2012h . 
Furthermore, the volatile content of habitable zone planets is 
expected to be highl y variable due to orbital migration (e.g., 
iRavmond et al.ll2004h . For instance, planets originating from the 
outer planetary system and made of a large fraction of cometary 
material can migrate to habi table orbital distan ces, resulting in 
the so-called ocean-planets dLeger et alj|2004t) . Planetary CO2 
reservoirs of the order of thousands of bars are certainly plausi- 
ble, when considering typical solar system values for the com- 
position of the cometary material. It is possible that the silicate 
carbonate cycle, which regulates the level of atmospheric CO2 
on Earth, does not operate on ocean planets in the absence of 
continents. Such large reservoirs of CO2 are therefore a concern 
for habitability if totally outgassed into a CC^-rich envelope. 

Fig.Q]shows the phase diagram of CO2. The critical point lies 
at r cr it=303 K and p c rit=73.8 bar. At a given surface temperature 
below r cl -i t , the vapor pressure curve actually limits the amount 
of CO2 which can be outgassed into the atmosphere, independent 
of the planetary reservoir. 
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Fig. 1. CO2 phase diagram. 

It is the aim of this study to quantify this maximum CO2 
partial pressure for a range of possible planetary scenarios near 
the outer edge of the HZ, based on the phase diagram in Fig. 
Q] In order to put constraints on atmospheric CO2, the inter- 
play between CO2 greenhouse effect, surface temperature and 
CO2 partial pressure must be taken into account. Therefore, this 
work will use an atmospheric model which consistently calcu- 
lates temperature profiles and surface conditions. It will be in- 
vestigated how different parameters such as planetary gravity, 
orbital distance and central star type, N2 pressure and surface 
albedo influence the behavior of the maximum CO2 partial pres- 
sure. 

The paper is organized as follows: Sect.|2]presents the model 
and scenarios, Sect. [3] the results and Sect. [4] a discussion. We 
conclude with Sect. 



2. Methods 

2.1. Atmospheric model 

We used a cloud-free, one-dimensional radiative-convective 
model to determine the globally averaged atmospheric 
temperature-pressure profile. 

The origin al model was first de scribed by IKasting et al.l 
d!984al) and IKasting et al 



were introduced by e.g. 



( 1984bl). Fur th er devel o pment s 



Kastind (Il988|). iKastingl dl99lh 



Kasting et al] dl993l)lMischna et al.N2000t) . IPavlov et all d2000b 



and ISegura et alj (2003b. The mode l ver sion used in this work 
is taken from lvon Paris et al.l (I2008I) and lvon Paris et al.l (l2010l) 
where more details on the model are given. 

The model atmospheres are assumed to be composed of N2, 
H2O, and CO2. Temperature profiles are obtained on 52 model 
layers, approximately spaced equidistantly in log (pressure). The 
pressure grid is determined from the surface pressure p sm f (vari- 
able, see below) up to a pressure of 6.6xl0~ 5 bar (fixed) at the 
model lid. 

The model calculates the temperature profile by solving 
the radiative transfer equation. The radiative fluxes are calcu- 
lated separately for the stellar (mostly visible) and the planetary 
(mostly IR) flux. The s tellar part of the ra diative transfer uses 
gaseous opacities from IPavlov et al] (l2000l) an d Rayleigh scat- 
tering formulations from lvon Paris et al.l d2010|). Gaseous opaci 



ties in the IR are based on Hitemp data dRofhman et al] 



and continuum absorptio n adapted from Clou gh et al 



1995) 



1989) 



and IKasting et al. (1984a). The purpose of the ID model used 
here is to calculate an arbitrary range of temperature-pressure 
scenarios ranging from the outer to the inner boundary of the 
HZ. Therefore, we used Hitemp in order to have reliable re- 
sults for wet, hot atmospheres. The choice of the specific opac- 
ity database (e.g., Hitran 2008, Hitran 2004, etc.) for gaseous 
absorption is not critical for the results presented below, i.e. rel- 
atively dry, cold scenarios. 

If the calculated radiative lapse rate is sub-adiabatic, the 
model performs convective adjustment, assuming a wet adia- 
batic lapse rate. This wet adiabatic lapse rate is determined con- 
sidering either CO2 or H2O as condensing species. 

The treatment of CO2 co ndensation for the calc ulation of the 
adiabatic lapse rate follows Ivon Paris et a i] d2010h . We assume 
that CO2 condensation occurs when the atmosphere is super- 
saturated with respect to CO2, as described by the super satu- 
ration ratio 5s: 



Pco, 

Pvap,C0 2 



= <S S = 1.34 



(2) 



where pco 2 is the partial CO2 pressure and p V ap,co 2 the 
saturation vapor pressure of CO2. The chosen value of «S S i s 
motivated by measurements reported in iGlandorf et al.l d2002l) . 
Condensation of an atmospheric constituent can occur when «S S 
is closer to unity than the value chosen here. N ote that other 
studies (e.g.. lKastingll99i1 or lKasting et al.ll993l) assumed «S S =1 
which represents the thermodynamic lower limit where conden- 
sation could occur. 

The water profile in the mo del is calculated based on th e 
relative humidity distribution of iManabe & Wetheraldl d 1967b . 
Above the cold trap, the water profile is set to an isoprofile taken 
from the cold trap value. Despite the fact that CO2 is allowed to 
condense, the major atmospheric constituents N2 and CO2 are 
isoprofiles throughout the entire atmosphere, i.e. are assumed 
to be well-mixed. The impact of fixing the CO2 mixing ratio 
at the saturation value on the atmospheric energy budget is ex- 
pected to be rather small, hence would not change our results 
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by much. A more consistent treatment of CO2 condensation (in- 
cluding an altitude-dependent CO2 profile) would involve verti- 
cal mass transport and an atmospheric pressure grid which is not 
in hydrostatic equilibrium in the region of CO2 condensation. 
Introducing this into our atmospheric model is beyond the scope 
of the current work. 

2.2. Model procedure 

The simulations started with a CO2 partial pressure of 73.8 bar, 
corresponding to the pressure at the critical point, p a -i t , and an 
isothermal temperature profile of 320 K, i.e. higher than the crit- 
ical temperature of 303 K. The choice of the initial temperature 
profile is not critical for the final outcome of the simulations. We 
did not allow for CO2 partial pressures higher than 73.8 bar, even 
though higher pressures are certainly possible (e.g., Venus). 

The surface pressure in model iteration step t + 1 is re- 
calculated based on the surface temperature r slll f as 

/>surf(Tsurf) = PN 2 + Pn 2 o(T sm i) + pCOafTsurf) (3) 

where is the fixed background pressure of N2. The water 
vapor pressure is obtained from 



ITERATION LOOP 



/>H 2 0(r s urf) - min(p vapj H 2 o(T sur f),Pocean) 



(4) 



with p V ap,H,o(^surf ) the water vapor saturation pressure at sur- 
face temperature and p ocea n the ocean reservoir assumed (here, 1 
Earth ocean, i.e. 270 bar). The CO2 partial pressure is accord- 
ingly calculated as 



Pco 2 ( T surf) = min(p vap ,co 2 (T sur f),pco 2 ) 



(5) 



Note that this corresponds to assuming a super-saturation ra- 
tio of S s -l at the surface, in contrast to *S S =1.34 used for the 
atmospheric CO2 adiabatic lapse rate (see eq. 0. This is mo- 
tivated by the fact that atmospheric condensation generally re- 
quires S S >1 (i.e., the presence of condensation nuclei). At the 
surface, however, atmosphere and reservoir are in equilibrium, 
hence the partial pressure follows the vapor pressure curve. 

The mixing ratio of N2 is then adjusted via 



C^ 



- C\ 



Psurf.t 
Psurf,t+1 



(6) 



where Cjv,«+i, C^ t are the N2 concentrations and /? S urf,t+i5 
Psurf.t the surface pressures at iteration steps (t + 1) and t. 

Based on the new value for the surface pressure p SUY f, the 
pressure grid on the 52 model levels is then re-calculated. 

Fig.|2]shows a flow chart of the model to illustrate the model 
procedure. 

The CO2 sa turation vapor pressure p V ap,co 2 is taken from 
lAmbrosel(ll956l) . It is divided into two temperature regimes. For 
T > 216.6 K (gas over liquid): 



fl?in(p va p,co 2 ) 
d\n(T) 



2.303 • T ■ 
( 867.2124 



(7) 



7-2 



+ 18.65612- 10" J - 



2 ■ 72.4882 • 10~° ■ : 
+3 ■ 93 ■ 10~ 9 r 2 ) 

For T < 216.6 K (gas over solid): 



isoprofile T=320 K 
Pco2=73.8 bar — 



Radiative 
Transfer 

Flux Profiles 



Radiative T profile 
Trad . 



convective adjustment 

from surface up to the 
tropopause 



Re-calculation of 
pressure grid 



P„,(f,.j\ baPW 



H ; profile: 

below cold trap: 



above cold trap: 
isoprofile 



Rel. humidity: R=R(p) 



Fig. 2. Flow chart of the model. 
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d\n(T) 



2.303 ■ T ■ 

1284.07 
(r-4.718) 2 



(8) 
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If surface temperatures remain above 303 K throughout the 
entire simulation, the maximum CO2 partial pressure is assumed 
to lie above the critical pressure. However, if surface tempera- 
tures converge to values below 303 K, the corresponding CO2 
partial pressure is taken as the maximum possible CO2 pressure 
for the particular planetary scenario. 

2.3. Parameter variations 

We varied five important model parameters: The planetary grav- 
ity, related to its mass and radius, the type of the central star and 
the energy input from the star, related to orbital distance, as well 
as model surface albedo and N2 partial pressure. Table Q] sum- 
marizes the varied parameters. 

- We assumed three different values for planetary gravity 
(lx, 2x and 3x Earth's gravity) which roughly corresponds 
to planetary masses of 1, 5 and 11 Earth masses, respec- 
tively, ac cording to mass-r adius relationships for rocky plan- 
ets (e.g.. lSotinetaT1l2007l) . 

- We used spectra of AD Leo, e Eri, the Sun and <x Boo 
as examples for M-, K-, G- and F-type stars, respectively. 
The same sample of stars has been used for numerous stud- 
ies regarding the influence of stellar t ype on atmospheric 



conditions (e.g., [S egura et al.l 120031 iSegu ra et 



mosphenc 
"all I20051 



Grenfell et a!ll2007al iGrenfell et alJl2007bl iKitzmann et al.l 
20101) . Stellar effective temperatures increased from M- to 
F-type stars, from about 3,400 K to 6,700 K, respectively. 
A more detailed description of the stellar s pectra as well as 
data so urces and references can be found in IKitzmann et al.l 
(120101) . 

The incoming stellar insolation 5/ at the top of the model 
atmospheres is calculated from 



S -S, 



(9) 



where So is the flux currently received by modern Earth (i.e., 
5'o=1366 WirT 2 ) and S is a constant factor related to orbital 
distance (e.g., for Earth, S = l). In this study, S was var- 
ied from 5=0.2 to 5=0.5. Corresponding orbital distances 
ranged from 0.2 1-0. 34 AU, 0.85-1.35 AU, 1.41-2.23 AU and 
2.67-4.22 AU fo r the M-, K-, G- and F-type stars, respec- 
tively (based on IKitzmann et al.ll2010l) . The range of stellar 
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Table 1. Parameter (range) for the runs performed 



Runs Gravity [g e ] Stellar type Stellar insolation S surface albedo As Pn 2 [bar] 



nominal 1-3 M-F 

N 2 study 1 M-F 

A s study 1 M-F 



0.2-0.5 0.13 1.0 

0.2-0.5 0.13 0.1-10. 

0.2-0.4 0.13-0.4 1.0 



insolation considered here roughly covers the outer limit of 
the HZ for the stellar types used in this work (e.g., GL 581 d 
with 5=0.29 and early Mars with 5=0.32, are both poten- 
tially habitable) as well as orbits slightly closer to or slightly 
farther away from the central star. 

- The above runs (nominal runs in Table [TJ were performed 
with the Nt partial pressure fixed at 1 bar. Increasing the 
amount of N2 (at fixed values of CO2 partial pressures) leads 
to two competing effects, a cooling effect (related to en- 
hanced Rayleigh scattering), and a warming effect (due to 
pressure broadening of absorption lines and continuum ab- 
sorption). Several studies have shown that increasing N2 par- 
tial pressures might indeed help to obtain h abitable surface 
conditions in atmospheric simulations (e.g.. iGoldblatt et all 
I2009L I von Paris et al 1 l2010h . Hence, we varied the N2 partial 
pressure from 0. 1 to lObar, for the Ig runs (N2 study in Table 

ED. 

- For all the model scenarios described above, the mea- 
sured rnean_sjrrfac^_jlbedoof the Earth (A sul f=0.13, taken 
from lRossow & Schifferl 1999b is used. However, the surface 
albedo has an im portant impact on th e calculated surface 
temperature (e.g., von Paris et al. l l2^lRosirigetini2010l 
IWordsworth et al] boiObTT Our model calculations do not 
take into account the possible increase of surface albedo due 
to condensing and freezing CO2 during the iterations. In this 
regard, our calculated CO2 partial pressures are likely to be 
upper limits. Measurements and m odeling of the albedo of 
CO2 snow by IWarren et alj {1990) suggest that the albedo 
of CO2 snow and ice might be significantly higher than 0.13. 
Therefore, we performed additional calculations (As study in 
Table[TJ with a surface albedo of A sur f =0.4 for the 1 g scenar- 
ios, at stellar insolations corresponding to 5 =0.2 and 5 =0.4 
and a N2 partial pressure of 1 bar, respectively. 



for the F- and G-star planets, CO2 partial pressures are below 
the critical pressure, at 10.9 and 23.2 bar, respectively. The upper 
stratosphere is sensitive to absorption of stellar radiation in the 
near-IR bands of CO2 and H2O, resulting in about 30 K increase 
for an M-star planet compared to the F-star planet. Additionally 
shown in Fig.|4]are the CO2 vapor pressure curve (<S S =1, dashed 
line, eq. |5j which intersects the temperature profile (for the F 
star and the G star) at the surface. Furthermore, Fig. |4] shows 
the CO2 condensation curve from eq.|2](>S s =1.34) indicating the 
CO2 convective regime. It is clearly seen that the atmospheres 
of the F-star and the G-star planet are dominated by a CO2 con- 
vective regime, followed by a very shallow near-surface H2O 
convective regime. In contrast, the K- and M-star planets show 
a relatively extensive lower troposphere dominated by H2O con- 
densation. 



S=0.35, 1g 




150 200 250 300 350 400 450 

T[K] 



3. Results 

Fig.[3]shows the maximum partial pressures of CO2 as a function 
of stellar insolation (hence, orbital distance, see Eq. [9} for the 
nominal runs of TableQ] Additionally shown as triple dot-dashed 
line in Fig.[3]is the CO2 partial pressure when using an equilib- 
rium temperature assuming zero albedo (i.e., the maximum equi- 
librium temperature, r eq inax , see eq.[TJ. This shows that detailed 
atmospheric modeling (taking into account the greenhouse ef- 
fect) is indeed needed to obtain consistent constraints on the CO2 
partial pressure. Also indicated in Fig.[3](by the horizontal plain 
line) is the boundary between liquid and solid phase of surface 
CO2, i.e. the triple point pressure of 5.1 bar (see the phase dia- 
gram, Fig. [TJ. For maximum CO2 pressures below 5.1 bar, the 
atmosphere is in equilibrium with CO2 ice, above 5.1 bar, the 
formation of (shallow) CO2 oceans is suggested. 

Fig.|4]shows sample temperatures profile of the simulations, 
i.e. a lg planet at 5=0.35, with a N2 pressure of 1 bar and 
As =0.1 3. As can be clearly seen, the K- and M-star planets re- 
tain their initial CO2 inventory of 73.8 bar (since at the surface, 
the atmosphere is not saturated with respect to CO2), whereas 



Fig. 4. Temperature profile for 1 g planets at 5 =0.35. 



3.1. Effect of stellar type 

From Fig. [3j it is clear that with increasing stellar effective tem- 
perature (changing stellar type from M to F), the maximum par- 
tial pressure of CO2 decreases. Also, the minimum stellar inso- 
lation 5 m in for which maximum CO2 pressures above p cr i t are 
possible depends sensitively on the stellar type (5 m j n =0.25 for 
the M-star planets and 5 m j n >0.5 for the F-star planets). This 
is due to the distribution of the stellar energy received by the 
model planets. With increasing stellar effective temperature, the 
stellar spectrum is shifted towards lower (bluer) wavelengths, as 
illustrated by Fig. [5] Broadly, the stellar spectrum can be sepa- 
rated into three regimes, 1) a Rayleigh scattering regime, 2) an 
absorption regime and 3) a "window" in between. The Rayleigh 
scattering regime is here defined as the spectral range where the 
Rayleigh cross section remains larger than 10~ 2 of the maximum 
value (A <0.75//m). 
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Fig. 3. Maximum CO2 partial pressure as a function of gravity and stellar insolation (as defined by Eq. [9j). The critical pressure 
Pcrit is indicated by the dot-dashed horizontal line. The triple-dot dashed line indicates the highest CO2 pressure calculated for the 
maximum equilibrium temperature (A=0, eq.Q]). 

partial pressures considered in this work, both the Rayleigh scat- 
tering regime and the absorption regime are almost entirely opti- 
cally thick to incoming stellar radiation (i.e., no radiation reach- 
ing the surface). In the Rayleigh scattering regime, radiation is 
reflected back to space (high spectral albedo), whereas in the ab- 
sorption regime, the radiation is deposited in the upper to middle 
atmosphere (very low spectral albedo), as illustrated in Fig.|6]for 
a 2 and 20 bar CO2 atmosphere. Depending on spectral type, the 
actual percentage of stellar radiation contained in the "window" 
changes quite considerably, as illustrated in Fig. [5] (around 50% 
for the M star, only 30% for the F star). 

Therefore, the planetary albedo becomes larger for increas- 
ing stellar effective temperature (M to F) because of the increas- 
ingly important contribution of Rayleigh scattering, and thus sur- 
face temperatures and corresponding CO2 partial pressures are 
lower. 




12 3 4 

wavelength [urn] 



Fig. 5. Cumulative energy of different central stars. Regimes are 
indicated by vertical lines. 



The absorption regime starts at about 1 .5 pm where the first 
strong water and CO2 absorption bands occur. At the high CO2 



3.2. Effect of planetary gravity 

The most noticeable effect when changing the planetary grav- 
ity g is the effect on atmospheric column density C. At constant 
pressure p, C and g are related linearly via C ~ pg~ l . Hence, an 
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wavelength [urn] 



Fig. 6. Spectral albedo for a 2 and a 20 bar CO2 atmosphere with 
surface temperature 288 K (corresponding to 17mbar of H2O), 
1 bar of N2 and As =0. 13. 1 g and 3 g planets indicated in black 
and red, respectively. Window regime is indicated by vertical 
lines. 



increase in gravity leads to a corresponding decrease of atmo- 
spheric column density. 

This leads to three important effects. Firstly, such a decrease 
in atmospheric column density leads to decreased Rayleigh scat- 
tering, hence a lower planetary albedo (see Fig. |6j, hence favors 
surface warming. Furthermore, less atmospheric column den- 
sity leads to less near-IR absorption of stellar radiation, hence 
higher albedo (again, see Fig. |6), hence surface warming (more 
starlight reaches the surface) and stratospheric cooling. On the 
other hand, a decreased atmospheric column density leads to less 
greenhouse effect (GHE), hence surface cooling. The net result 
on surface temperature when combining these three effects (ei- 
ther cooling or warming) depends on the amount of CO2 and the 
stellar type which determines the planetary albedo and stellar 
energy distribution (see Figs.[5]and[6]l. 



Fig. [7] shows the ratio between calculated CO2 pressures at 
Ig and 3g. At low stellar insolation, hence low CO2 pressures 
(see Fig. [3}, increasing gravity leads to cooler surface tempera- 
tures, and consequently lower CO2 partial pressures (i.e., a ra- 
tio higher than 1 for all stars except the F star in Fig. |7). This 
indicates that the impact of the reduced GHE is dominating, in 
agreement with other studies of opt ically rather thin planetary at- 
mospheres (e.g.. lRauer et al.luOllh . In contrast, at higher stellar 
insolation (and correspondingly higher CO2 pressures), increas- 
ing gravity leads to warmer surface temperatures, hence higher 
CO2 partial pressures (i.e., a ratio lower than 1 in Fig. [7]), imply- 
ing that the decrease of the GHE is compensated by the decrease 
in planetary albedo. The influence of the stellar type is clearly 
seen in Fig. [7] For the M-star planet, with very little radiation in 
the Rayleigh regime (see Fig. [5]), the effect of increasing grav- 
ity is much higher than for the F-star planet, for which Rayleigh 
scattering is very important. 



3.3. Implications for habitability 

As can be inferred from Fig. [3] our calculations imply that rela- 
tively massive CO2 atmospheres of the order of several bars are 
possible for almost all scenarios, even for planets orbiting far 
from their central star (stellar insolation S >0.25). 

At the triple point temperature of water, i.e. 273 K, which 
permits its liquid phase, the CO2 vapor pressure is about 34 bar 
(see Fig.Q]). Hence, Fig. [3] implies that liquid surface water can 
be achieved for stellar insolation S 34^- as low as S34bar— 0.25 for 
the M-type star and 534bar=0.4 for the F-type star, providing a 
sufficiently large source of CO2 is available for outgassing on 
the planet. This is, however, not the outer edge of the HZ, since 
surface temperature is not necessarily a monotonic function of 
CO2 part ial pressure (known as the maximum greenhouse ef- 
fect, e.g jKasting et al.ll993l) . The CO2 pressures corresponding 
to the maximum surface temperatures are therefore expected to 
be somewhat lower than the maximum CO2 pressures in Fig. [3] 
Hence, the outer edge of the HZ is most likely located at lower 
stellar insolation (i.e., farther away from the star), than S^bar- 
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Fig. 7. Effect of varying planetary gravity on the calculated max- 
imum CO2 pressures. Ratio between calculated CO2 pressures at 
lg and 3g. (Super-)critical pressures which have a ratio of unity 
at higher values of S not shown. 



3.4. N 2 partial pressure 

The results of the N2 study (Sect. I2.3I and Table [TJ are shown 
in Fig. [8] As expected, for the high CO2 partial pressures found 
for higher stellar insolation, the effect of varying N2 is negli- 
gible, given that CO2 is a much more efficient Rayleigh scat- 
terer than N2. However, for lower stellar insolation, and conse- 
quently lower CO2 partial pressures, the effect of N2 becomes 
discernible. 

At these lower stellar insolation, the warming effect of 
adding N2 to the atmosphere is clearly dominating, since the 
calculated maximum CO2 pressures increase with increasing N2 
partial pressure. The effect is rather pronounced (almost a fac- 
tor of 4 when increasing from 0.1 to lObar) for the M star 
since Rayleigh scattering does not contribute greatly to the over- 
all energy budget for these cases (most of the stellar radiation is 
emitted at wavelengths where Rayleigh scattering is negligible, 
see Fig . [5j ■ F° r the F-star simulations, maximum CO2 pressures 
increase only by about 30%, i.e. warming and cooling effects 
approximately cancel out. 
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Fig. 8. Effect of varying N2 partial pressure p^ 2 on the calculated 
maximum CO2 pressures. (Super-)critical pressures which have 
a ratio of unity at higher values of 5 not shown. 



3.5. Surface albedo 

The results of the surface albedo study (Sect. I2.3I and Table [TJ 
are presented in Fig. [9] which shows the decrease in calculated 
maximum CO2 pressure when increasing the surface albedo. At 
5=0.2, the decrease of CO2 pressure is rather large, reaching 
about a factor of 20 for the M-type star. For a planet orbiting 
around an F-star, calculations imply maximum CO2 pressures of 
the order of 0. 1 bar, so rather a teneous atmosphere. At 5 =0.4, 
the effect of increasing surface albedo is smaller than at 5 =0.2, 
but still reaches about a factor of 2-3 for the F-type star. 
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Effect of varying surface albedo on the calculated max- 
CO2 pressures. (Super-)critical pressures at 5=0.4 for 
and K star not shown (indicated by horizontal line at 
73.8bar). 



Fig. [9] shows that, at 5=0.2, the M-star planet is much more 
sensitive to a change in surface albedo (a reduction of a factor of 
about 20 in CO2 pressure) than the F-star planet (a factor of 8), 
as seen by the steeper slope of the M-star line. The sensitivity is 
generally increasing for increasing stellar effective temperature 
(type from M to F). This is due to the larger amount of stel- 
lar energy emitted in the window regime (see Sect. 13. H and Fig. 
[3}. Hence, the response to an increase in surface albedo, which 
affects principally the window, is more pronounced for the M- 
star planet and for lower stellar insolation (and correspondingly 
lower CO2 partial pressures). For example, at 5=0.4, the reduc- 
tion for the F-star planet is decreased to about a factor of 2. 

In order to investigate the combined effect of, e.g., an in- 
crease in N2 partial pressure and an increase in surface albedo, 
we performed some additional test runs with both parameters 
changed. For the M star case, for example, the effect of N2 was 
nearly unaltered even at high surface albedo. At 5=0.2, an in- 
crease in surface albedo reduced the maximum CO2 pressure 
from 8.2 to roughly 0.4 bar (see Fig. |9]l whereas an increase 
of N2 partial pressure increased the maximum CO2 from 8.2 
to 16.3 bar (see Fig. [|J. At high surface albedo and high N2 
pressure, the maximum CO2 pressure obtained was 14.3 bar, i.e. 
nearly as high as for the simulations at low surface albedo. 

4. Discussion 

4.1. H2O-CO2 oceans 

As has been shown above (Fig. [3}, for planets orbiting within 
the HZ of K-G-F stars there is a region of liquid surface CO2 
combined with surface temperatures above 273 K, i.e. liquid 
surface H2O. This means that it is possible to form H2O- 
CO2 oceans. Then, the question of planetary habitability would 
depend strongly on the pH of the liquid, even though ex- 
tremophiles on Earth could su pport quite low pH values (e.g., 
iRothschild & Mancinellill200lh . A detailed investigation of this 
interesting issue is however beyond the scope of this work. 

4.2. Implications of model assumptions 

The ID atmospheric model used in this work is based on rela- 
tively few, simple assumptions. Most of these assumptions are 
physically justified, i.e. the assumption of adiabatic temperature 
gradients in the troposphere or radiative transfer as the main en- 
ergy transport mechanism in the upper atmosphere. However, 
some of them (presence of clouds, greenhouse gases, water pro- 
file, etc.) are model-specific, hence need to be discussed further 
with respect to their possible influence on the results presented 
above. 

The model is a cloud-free code, hence the poten- 
tial impact of CO2 clouds on the climate is neglected. 
It was shown by several au thors that this potentia l im- 



pact could be quite large (e.g. , Forget & Pierrehumb ertl 1 1 99 



Mischna et aT1l2000i IColaprete & Toonll2003l IWordsworth et al 



2010bi rWordsworth et alj|201 ll) . However, this effect depends 



sensitively on cloud opacity, cloud coverage and cloud altitude. 
In addition, the effect of clouds is also probably very dependent 
on stellar type (see, e.g. jKitzmann et ai]|2010l investigating the 
effect of stellar type for H2O clouds). Investigating this is there- 
fore a subject of further studies. 

Furthermore, the model atmospheres considered in this work 
contained only the greenhouse gases CO2 and water. This choice 
may be restrictive when applied to our own solar system, since 
other species, such as O3, SO2, CH4, and N2O, have been con- 
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sider ed in models of t he early Ear t h or early Mars climat e 
(e.g., lYung etall 119971 iBuicl 120071 iHaaa-Misra et afl l2008h . 
But given that the concentration of these gases depend on very 
specific planetary scenarios (e.g., outgassing history, biospheric 
evolution, etc.), assuming them in the context of exoplanets 
(without any geological or other constraints) is rather arbitrary. 
However, the impact on stratospheric temperatures through the 
absorption of UV (e.g., O3 and SO2) or near-IR (e.g., CH4) stel- 
lar radiation is potentially important. 

Radiative transfer in dense, CC^-dominated atmospheres 
presents many challenges (e.g., collision-induced absorption, 
sub-Lorentzian behavior of line wings, etc.). The parametriza- 
tion of the collision- induced absorption (C IA) used in this 
study is taken from Kasting et al.l (Il984bl) . A recent study 
dWordsworth et alj|2010al) presented a re vised parametrization , 
showing that the calculation presented bv lKasting et al.l d!984bl) 
most likely over-estimates the opacity. In order to estimate the 
impact of the CIA uncertainties on our results, we performed 
a sensitivity study with a reduced (by roughly a factor of 2) 
CIA. The conclusions however did not change qualitatively. At 
5=0.2, calculated CO2 maximum pressures around K-, G- and 
F-stars decreased by less than 50%, for the M-star the maximum 
CO2 pressure decreased from 8.2 to 3.0bar. At S =0.35, results 
changed less than 20% except for the K-star planet, where a 
maximum CO2 pressure of 47.3 bar was calculated, instead of 
73.8 bar (i.e., the critical pressure of CO? , see Fig.[3]>. Therefore, 
our calculations (using Kas ting et alj 1 1984bh are likely to be 
overestimates of the maximum CO2 partial pressures. 

The model uses a super-saturation of <S S =1.34 to determine 
the CO2 convective regime (see eq.[2j. The choice of S s has been 
sho wn to be very impor tant for early Mars climate simulations, 
e.g. lPollack et al.l d 1987b (using S s =<x>) find signi ficantly higher 
surface temperatures (>30 K) than Kas ting! (1 1 99 ll) (usin g <S S =1). 
The assumed £,=1.34 is based on lGlandorf et a TT d20Q2h . a value 
observed for specific conditions (e.g., dust loading available for 
nucleation) which could be different on exoplanets (as low as 
Jvs=l, but also possibly significantly higher). In this sense, the 
calculated maximum CO2 pressures are not necessarily upper 
limits. To further investigate this, we performed some sensitiv- 
ity simulations with S s =l. As expected, calculated maximum 
CO2 pressures were lower, of the same order of magnitude as 
for the CIA study mentioned above. However, the main conclu- 
sions obtained in this work (i.e., the existence of maximum CO2 
pressures far below the critical pressure) were not affected. 

The relative humidity profile used in this work 
(iManabe & Wetheraldl 1 19671) has been derived from observa- 
tions of modern Earth. It has been used in many I D simulations 
of terrestrial exoplanets, both Earth-lik e (e.g.. | Segura et"al] 
20031 iGrenfell et al.ll2007al) and not (e.g.. lvon Paris et al]l2010l 
Wordsworth et al.l2010bl) T Since the humidity profile is anything 
but trivial to model in ID simulations, some authors c hose to fix 
relative humidity at an isoprofile (e.g.. lKasting|ll99il) . However, 
given the large amounts of CO2 in the model atmospheres 
(73.8 bar at 303 K), the impact of water (42mbar at 303 K) 
on atmospheric structure (via near-IR absorption) and surface 
conditions (via the GHE) is somewhat negligible. Therefore, the 
choice of the relative humidity profile is probably not important. 



4.3. Synchronous rotation 

For planets orbiting very close to their star, tidal locking of 
the planetary rotation with the orbital period is very likely. The 
time scale fi oc k of tidal locking is very sensiti ve to orbital dis- 
tance (fi oc k ~ a 6 , a orbital distance, see e.g. iGriefimeier et alJ 



120051) . Hence, tidal locking is mainly an issue for the habit- 
ability of planets orbiting around M stars due to the closeness 
of the HZ to the star. It has been argued that for planets with 
a perpetual nightside, the atmosphere could collapse since the 
nightside forms a cold trap for the volatiles, which, in the con- 
text of this work, could present an alternative way of obtain- 
ing maximum CO2 pressures. Ho wever, as has been shown by 
numerous modeling studies (e.g.. | Joshi et al.lll997l I Joshil [2003L 
IWordsworth et al.ll2.01 ll Kite et al. 201 1 ). moderately dense at- 
mospheres containing hundreds of millibars or more of CO2 
are sufficient to avoid atmospheric collapse by means of atmo- 
spheric circulation. Hence, the M-star simulations presented in 
this work are not thought to be subject to atmospheric collapse 
induced by synchronous rotation. 

5. Conclusions 

We have presented a detailed parameter study to constrain the 
maximum CO2 partial pressure possible for terrestrial exoplan- 
ets, using a ID cloud-free atmospheric model. Parameters inves- 
tigated included the central star type, the orbital distance and the 
planetary gravity. Furthermore, we investigated the influence of 
N2 partial pressure and the surface albedo on the maximum CO2 
partial pressure. 

Results imply that super-critical atmospheres (i.e., 
Pco 2 ^Pcrit =73.8 bar) are possible for planets around M 
stars for stellar insolation corresponding to 5^^=0.25 or higher. 
For increasingly bluer stars (i.e., higher effective temperatures), 
this super-critical stellar insolation increases (e.g., 5 , cl i t >0.5 for 
an F-type star). For lower stellar insolation, the calculations 
presented here imply that there is indeed a maximum CO2 
partial pressure, even if the planets are orbiting well within the 
habitable zone. Nevertheless, massive CO2 atmospheres of the 
order of bars are still possible for most scenarios. For planets 
orbiting very far from an F-type central star (e.g., 5=0.2 in this 
work), CO2 partial pressures could be constrained to be less 
than 1 bar. 

The effect of planetary gravity is twofold. For low stellar in- 
solation and corresponding cold surface temperatures, increas- 
ing planetary gravity leads to a decrease of maximum CO2 
partial pressure due to less atmospheric greenhouse effect. At 
higher stellar insolation, an increase of planetary gravity in- 
creases the calculated maximum CO2 partial pressure because 
of less Rayleigh scattering. For F-star planets, the effect is up 
to a factor of 2, depending on stellar insolation, whereas for an 
M-star planet, a factor of about 4 has been calculated. 

Increasing the N2 partial pressure leads to warmer surface 
temperatures for all cases, hence corresponding maximum CO2 
pressures are higher. The effect reaches up to a factor of 3 for 
planets around an M star, upon increasing the N2 partial pressure 
from 0.1 to 10 bar. The surface albedo has an important effect 
on the values of the maximum CO2 partial pressure. A higher 
surface albedo leads to cooler surface temperatures, hence less 
CO2 in the atmosphere. Decreases of about a factor of 20 have 
been shown when increasing the surface albedo from 0. 13 to 0.4. 

The presence of CO2 and H2O clouds could alter these re- 
sults because of their potentially large impact on the planetary 
energy balance. However, our (clear-sky) results show in a ro- 
bust way that the composition and the evolution of planetary 
atmospheres strongly depend on orbital and planetary parame- 
ters. Although a consistent model for the determination of CO2 
partial pressures must take processes such as sequestration or 
outgassing into account, our results show that there is a funda- 
mental thermodynamic limit to the amount of CO2 in terrestrial 
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atmospheres, independent of the planetary reservoir. Hence, a 
more detailed coupling between interior, surface and atmosphere 
models should be used to accurately predict atmospheric compo- 
sition of terrestrial planets. 
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